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Production of Dirac particle in a deformed Minkowsky space-time 
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In this paper we study the Dirac field theory interacting with external gravitation field, described 
with tetrad of the form e^ ( x ) = 4- co£ a x a ), where e = 1 for p = 0 and e = i for p = 1, 2,3. The 

probability density of the vacuum-vacuum pair creation is given. In particular case of vanishing 
electromagnetic fields, we point out how this deformation modify the amplitude transition. The 
corresponding Dirac equation is solved. 


I. INTRODUCTION 


The Dirac particle theory arise in the theoret¬ 
ical description of the fermion particles phenom¬ 
ena. They also become central to elementary par¬ 
ticle physics, as the starting point for quantum the¬ 
ory of electromagnetic interaction. In the past few 
years this theory have been widely studied in various 
curved backgrounds due to its importance in both 
astrophysics and cosmology, as well as in the study 
of particle creation processes [l]-[6] . Pair produc¬ 
tion is a phenomenon of nature where energy is con¬ 
verted to mass, nevertheless there are only few prob¬ 
lems for which the Dirac equation can be solved ex¬ 
actly. Some of them are give in [7]-[15] and references 
therein. The explicit solution is crucial in the parti¬ 
cle creation processes and is the base of the standard 
cosmological model [25]-[27] . This amounts to claim 
that, the formulation and behaviour of fermion par¬ 
ticles physics including the gravitation field is per¬ 
formed using solution of Dirac equation. 

In the present work, the transition amplitude and 
the probability density of the pair creation of the 
Dirac particle is examined in the twisted Minkowsky 
space-time. A solution of the Dirac equation is pro¬ 
posed. The metric is chosen to be the first-order 
fluctuation of the flat Minkowsky pseudo-metric ten¬ 
sor rj ab . We pointed out how this deformation of flat 
metric modify the well know probability density of 
creation of Dirac particle, a while ago computed in 
[2], [4] and references therein. 

The paper is organized as follows. In section (II), 
we quickly review Dirac particle theory, interacting 
with gravity field. In (III) we compute the probabil¬ 
ity density of pair production. The case of vanishing 
electromagnetic (EM) fields is examined. The solu¬ 


tion of the corresponding Dirac equation is proposed 
in the section (IV). In the last section (V), we con¬ 
clude our work and make some remarks. 


II. DIRAC EQUATION COUPLED WITH A 
WEAK GRAVITATION FIELD 


In the curve space-time the conventional affine 
connection is replaced by the spin connection 
T^ which is expressed in terms of the vierbein fields 
( e%(x )) (see [16]-[24]). Then any curved space de¬ 
scription of physics can be replaced by an equivalent 
and simpler flat space physics, through the vierbein 
transformation. So there was an equivalent formu¬ 
lation of general relativity involving the dynamics of 
the so-called spin-connection. This approach came 
to be known as Einstein-Cartan theory and leads to 
consider general relativity as gauge theory approach 
of gravity [18]. 

An arbitrary geometrical object defined on the 
Riemann space-time manifold can be locally pro¬ 
jected on the tangent Minkowski space, simply by 
contracting its curved indices with the vierbein and 
its inverse. For rank n tensor object T, we can write 


T a I 02 '" 0 " = e£ ( x)e% (*)■•• e£ , (1) 

T cna 2 -a n = e ai { x ) e al ( X ) ' ' ' e aZ , ( 2 ) 

where the Latin indices (a, b, c, ■ ■ ■) is used only 
for the flat space-time and the Greek indices 
(a, /3, p, ■ ■ ■) for the curve space-time. The “tilde 
notation” is used only for the curve space variables. 
e“(x) represent the inverse of e%(x). The metric ten¬ 
sors g M " and r] ab = diag(l, —1, —1, —1) are related by 
g^(x) = e%(x)e%(x)ri ab , or r] ab = e^(x)e b (x)g 111 ' (x). 
The connection T^ is 



where a 1 * 1 ' = 


^[7 li ( x ),'y v { x )]- 
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We consider the Dirac equation coupled with both 
gravitational and EM fields, given by the following 
relation 

(vy ll (x)D IJj - m)ip{x) = 0, (4) 

where D M + iA ll . In this expression, 

the vectors T^ and A^ are respectively the gravi¬ 
tation and EM gauge vectors. The field ip, is a four- 
components complex functions of space-time coordi¬ 
nates x M , p = 0,1, 2,3. The Dirac gamma matrices 
7 M (:r), acting on the vector fields ip, satisfy the anti 
commutation formula: 

{ 7 M (cc), 7 I/ (a:)} = 2g M " such that 7 M (®) = e%(x)j a .(5) 


such that the metric depend only on the coordinates 
(t = x°, x = x 1 ), i.e. 

uj 2 = uj 3 = UJq = uj, Ujf = uj. ( 10 ) 

The vector T M is 


/ r° \ 

f f iV \ 

Fl - 

f 7 °7 1 

t 2 - 

f 7°7 2 - f 7 V 

\r 3 J 

\ f'7°7 3 ^ / 


We get the following result: 

Proposition 1. The Dirac equation in the curve 
background defined with the metric (9) and coupled 
with EM fields A^ = (0, 0, Bx, — Et ) is given by 


The flat space-time gamma matrices 7 are expressed 
with the Pauli matrices a 1 , i = 1, 2,3 by 


7 


0 


1 2 0 \ i_f 0 a 1 
0 - 1 2 ) ’ 7 ^ -a 1 0 


1 , 2 ,3.(6) 


Remark that the equation (4) provided from the 
Euler-Lagrange equation of motion of the action S : 


S = 


d i xyf ZI g (iipY { x )Dnip 



( 7 ) 


where ip = ip^°(x). 

The dynamics described by the relation (4) is in¬ 
variant under an external local transformation of 
Lorentz group. In the case of internal local trans¬ 
formation this invariance is satisfy for the dynamics 
occurring within the space-time manifold. The sym¬ 
metries of this internal space are chosen to be the 
gauge symmetries of some gauge theory, so a unified 
theory would contain gravity together with the other 
observed fields. 

For \ujf\ << 1, we consider the vierbein field e%(x) 
as 


K)( x ) = dia § l+w°x a , 7(1 + 07 ( 2 ;“), 

i{l+u 2 a x a ),i{l+ulx a )\. ( 8 ) 


We shall use the notation cu%x a =: oj£ a x a . Then, the 
metric tensor g^" =: 77 ^" + (where is the 
perturbation tensor), takes the form 


(gH = diag[l + 2uj° a x a , -1 - 2 uj l a x a , 

-l-2u;;y,-1-2^1, (9) 


such that the limit where (w) —> 0 restore the 
Minkowsky pseudo-metric, g^ 1 " can be considered 
as the first-order fluctuation of the flat Minkowsky 
pseudo-metric. Now let us choose the tensor (w) 


ry°<9o — 7 ■’dj — i"f 2 Bx + i'f’Et + 9 * 077 ° — 3 W 7 1 
—m(l — cot — uix) ip(x fJ ‘) = 0, j = 1, 2, 3. (12) 


The solution of the corresponding equation can be 
split into 

ip(t, x, y, z) = ip(t, x) exp [i(k 2 y + k 3 z )\, (13) 

where ip(t,x) is a function which depends only on t 
and x. 


Proof. Using the relation ( 8 ), the Christoffel tensors 
are: 

Ka = -<*>“, r b ab = -UJ b a , 

^aa = VaaVbbOJb, Kb = 0, a ^ b ^ C, (14) 

where the Einstein summation are not taking into 
account in the above relations. Also the components 
of the Lorentz connection are 

To = \ (<A 01 + cA 02 + co°a 03 ), (15) 

Ti = -^ 0 V 10 -^ 12 -^ 13 ), (!6) 

r 2 = --[iu 0 a u - o+cr -o^cr^J, (17) 

ti _ 1 (■ ,3,30 ,,3,31 , ,3 3 2 \ /ini 

f 3 _ 2 V W ° ct - u 1 a — uj 2 a J, (18) 

which are reduced to ( 11 ) using ( 10 ), and a ab = 
d [ 7 a , 7 b ]. Now, by replacing the expressions (15) in 
(4), the Dirac equation becomes 

^7 M (x)5 /i -7 /i (x)A M +w a 7 a -m^V , (* M ) = 0, (19) 

where ujq = *(wq + uj 3 + w o) = 3*oi, uj i = (a;) 1 — uif — 
w i) = -w, uj 2 = (w^-w^-wf) = 0 , w 3 = (w§-w|- 
uj 3 ) = 0. We choose the external electromagnetic 
field as E = Ee x , B = Be x , where e^, is the unit 
vector in x direction. One solution of the Maxwell 
equation is then A^ = (0,0, Bx, —Et). Finally, the 
relation ( 12 ) is well satisfy. tj 
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III. TRANSITION AMPLITUDE OF THE 
MODEL 


+ (4 7 °7 3 - AEt 2 )P 3 - 67°7 3 Et + 4 7 °7 2 BX 
+ 2ij 1 'y 2 Bt + 2t(B 2 X 2 + E 2 t 2 ), (26) 


In this section we study, how this new met¬ 
ric modify the pair creation of fermion particles. 
We consider the Hilbert space of coordinates vec¬ 
tors "H x such that the space-time coordinates = 
(x°, x 1 , x 2 , x 3 ) = ( t,x,y,z ) are eigenvalue of coor¬ 
dinate operators = (H4, A, Y, Z) acting on P x , 
i.e. for | t,x,y,z >€ P x 

X^\t,x,y,z >= x^\t,x,y,z > . (20) 

The Hilbert space of momentum space vectors TL P 
is define as the Fourier transformation of TL X . The 
momentum operator P M = (Pq, Pi, P 2 , P 3 ), (Po = 
id 0 , Pj = —idj, j = 1, 2, 3), is defined by 

P?\ko, ki,k 2 ,k 3 >= k„\k 0 , ki,k 2 , k 3 >, 

gikpX 1 * 

< k 0 ,ki,k 2 ,k 3 \t,x,y,z >= N e K.(21) 


y 2 H = 2A|P| 2 - 27 1 7 0 P 0 + (2i7 1 7 2 + 4HA' 2 )P 2 
+ (2i 7 1 7 3 - 4PAf)P 3 - *7 1 7 3 Pt + 1 7 2 PA 

- 2 r y°"f 3 EX + 2X(B 2 X 2 -)- E 2 t 2 ). (27) 

Then V =: \Z(A,Y)\ 2 is explicitly written as 


V = exp 


— Tr In 


£(u, u>) 


= exp 


— Tr 


1t( 


n(u>, w, E, B). 


_. gisn(u),E,B) _ ^is£(uj 




(28) 


where 


£(oj,u) = X H {0,0) + X h (oj, 0) + X H (0,u), (29) 


n{u>, u>, E, B) = y H {0,0,E,B) +y H {u,0,E,B) 
+y H (P,u,E,B). (30) 


Also the curve space-time coordinates operators are 
given by A^ = ( t,X,Y,Z ) and conjugate momen¬ 
tum operators P M = (Po, Pi, P 2 , P 3 ) such that Po = 
*good 0 and Pj = *g jjdj, j = 1,2,3. 

In the path integral point of view, the action (7) 
gives the transition amplitude of the model (or the 
partition function Z(A,T) = A f J DifDijj e lS ) which 
is explicitly written as : 


Z(A, T) = exp 


Tr j ~m + ic 

M 


( 22 ) 


with M = *7 m (9 m — T^ + iAy) — m + ie and 
the normalization constant is defined such that 
21(0,0) = 1. Using the followings identities: 

ClaC- 1 = -7*,, (rY = -e^C- 1 = 

-C^C- 1 , T^ = —CT^C -1 , we come to M l = 
iC^^{d^ — T m + iA^C -1 — m + ie and then the 
conjugate of the functional Z(A,Y) is given by 


Z t (A, T) = exp 
where C = Ppy 0 . 


- Tr In 


iC^C 1 d^ l + m + ie 





We now compute the transition amplitude 
\Z(A, T)| 2 . For this, let us define the quan¬ 
tities Xh{oj, 0) = ojXfi, Xh{w, 0) = ul’g, 

y H (u,o,E,B) = uyj I , y H (o,u,E,B) = zy 2 H , such 

that 


X l H = 2t\V\ 2 + 7 °7 1 P 1 + 7°7 2 P 2 + 7 °7 3 ^3, (24) 


We get the following statement: 

Proposition 2. Consider that the EM fields are 
vanishing. For very small positif parameter e of the 
size 1/w 2 , the probability of the pair production takes 
the form 


V = exp 


7rMe^m 8 r r 37 

- 2 N~ - 

1024ve L 7 12 


1 

3 


(4 In b -In a) 


( 3 !) 

where u 1 = vu, a = |^(v 2 + 1), b = |^(4t) 2 + 1), 
A 7 is the Euler number, M = f y,(y, z)dy dz, and 
p{y, z) is the test function. 


Remark 1. Note that the case where M > 0 is not 
fulfils. This leads to a infinite probability density. 
We choose the test function y(y, z) such that M < 0, 
and then 


Tr MeWr 37 

_ 0 N _ 

1024i> l_ 7 12 


1 

3 


(4 In b — In a) 


> 0. 


In the figure (1) we give the plot of this probability 
density as function of the parameter e. This figure 
gives asymptotically the values of /(e) = V when e 
tends to zero. Then the limit e —> 0 leads to V ~ 0. 


Proof, of the proposition (2): The rest of this section 
is devoted to the proof of the proposition (2). The 
density Yl =: YI(oj,uj,E, B) such that V = exp(— Yl) 
can be expanded as 


H = Yl\{w,u),E,B) — Yt 2 (<jj,w,E,B), (32) 


X 2 = 2A|P| 2 + 7 °7 1 Po + tyVPa + (25) 

yjj = 2f|P| 2 + 47 0 7 1 P L + (47°7 2 +4PAf)P 2 


where 

Yh = j dx^°° ^(xle^^’^lx), 


(33) 


3 




















( 40 ) 



FIG. 1: Plot of V = /(e) with v = 1, M = —2n, m = 1, 
00 = \ft. 


= dx 


ds 


X |e is ^(“’ S )| X >. 


(34) 


We consider the mean values: 

Xh = <x|A^|x>, x 2 H =<x\Xjj\x>, 

Vh = < x lJ4l x >> Vh = < x l(yfrl x > • (35) 

Now we choose u> = vuj and vt + x > 0 and we are 
focussing on the computation of fii (w, oo, E, B). We 
get 


n, = fi 


■ x | e <«ya(0,0,B,B)| x v 


u 

j dxdk ^1 + isuoy+ isdoy^J (36) 


with 


x e 


^|x> = 


—iEB coth(.Es) cot(i3s) _ ism 2 

4^ 6 


A simple routine checking shows that 


(37) 


dxdk exp 


VI ne 
4 s 2 e 


( 1 

/ dxdt ( 
(vt + x) 2 l 


isujlvyn+y 2 H 


exp 


ISLU I — 


(vt + X) 2 


with 


(38) 


P(x,t) = 2u 2 (3 + 2Bxtj°"f 2 + B 2 x 2 t 2 

- 2j°'y 3 Et 2 + E 2 t A ) + 2v(2Bx 2r y 0r y 2 
+ i^rfBxt + 2B 2 x 3 t - 2'^'f‘Ext 

- i'ir'Et 2 + E 2 xt 3 ) + | + 2i"/ 1 'y 2 Bx 2 

+ 2B 2 x 4 - 2ij 1 j 3 Ext + 2E 2 x 2 t 2 , (39) 


and 


Q(x,t) = u[ 4 7 V^- 6 7 °7 3 ^ + 2n 1 7 2 Bt 


+ 2 t(B 2 x 2 + E 2 t 2 )\ - iry'-'fEt 
+ 3* 7 1 7 2 i?a: — 2 7 ° 7 3 Ex 
+ 2x(B 2 x 2 + E 2 t 2 ). 

In the same manner f lo takes the form 


<u= r dx *k 


JsX H (0,0) | 


■'ll " 

x (l + isuox^j + isoox 2 ^ 


where 


(x|e Is ^|x) = - 


167t 2 s 2 


(41) 

(42) 


We can now show that 


dxdk (1 + isvooxh, + isoox 


L H I 


= dx 


4 s 2 e(vt + x)‘ 


exp 


3 isuo(v 2 + 1) 
8 (vt + x) 


However 


Q(uj,E,B) = 


iMe 71 
64e 


ds ■ 


—4EB coth(£s) co t(Bs)J(to, Xo) 


(43) 


oX(*0, X 0 ) 

S z 

(44) 


where M is chosen to be M = / ^(x, y)dydz < oo, 


i(t Q ,xo) = r [ 

J Xn J ti 


dxdt 


7re 


Xq j to 

exp 


4 s 2 e(vt + x) 2 


3 isuj(v 2 + 1) 
8(i?£ + x) 


n 


(45) 


and 


/»oo /*oo 

J(t 0 ,x 0 ) = / 

J xn Jtc 


dxdt ■ 


1 


x exp 


so Jt 0 + x ) 2 

P(x,t) 

isuo | - —-r + <y(a;, t 

(vt + x) 


A 


(46) 


The integral (44) exhibit the divergence at point x = 
t = 0. This shall be regularized by using the Cauchy 
principal value. For E = B = 0 we get 


1(0,0) = 


dt 


a i Jo 

in 1 
2v v . 

pOO 


sm(ati/vt) + 2isin 2 (ai/2ut) 


1 — AL — In — 


(?) 


(47) 


dt 


sin(a 2 /uf) + 2isin 2 (a 2 /2uf) 


1 - AL - In — 


(48) 


J( 0,0) = — / 

a 2 Jo 
in 1 

2v v - 

with 07 = '^(v 2 + 1), a 2 = ^r(12u 2 + 3) and IV 7 is 
the Euler number given by Nj = 0.577215664. Also, 
for o/el the integral 


(?) 


Q = 


roo -ism 2 . . 

J ds - g — ^4 ln(6s) — ln(as) J 
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admits the Cauchy principal value 


/C2 satisfying the commutation relation [/Ci, AH2] = 0 
and given by 


Pv(Q) = — [415 - 3001V 7 + 72V 2 + 12tt 2 

+ 961n6 2 — 24 In a 2 — 3001n(im 2 ) 

+ 72 ln(im 2 ) 2 + 144iV 7 ln(im 2 ) 

+ 4^ — 25 + 12TV 7 + 121n(W)) In a 

- 16 ( - 25 + 12 V 7 + 12 ln(*TO 2 )) In b 

(49) 

Remark that the probability density of pair creation 
in the limit 25 = 0 (see [2]) correspond to 


= ff £ (coth(fa|)exp( 


“E 1 

fc=l 

Using the Taylor expansion as 


knur 

iT" 


(50) 


12(25, V, B) = 12(0, E, B) + 2512'(0, E, B ) + 0(25 2 ), 

(51) 

we come to 12 (w, 0, 0) = 2512'(0, 0,0) + 0(u 2 ) and 


7T poo „—ism 


iA/fp , 71 1 w ~ e 

12(25,0,0) = —^/ ds 6 —jr—T(0,0), (52) 

64e Jo s 5 

r(o,o) = j(o,o)-4j(o,o). 

We choose the real part of 12(25,0,0) denoted by 
3? e 12(25,0,0) = [12(25,0,0) + 12* (25,0,0)]/2. Using 
(47), (48) and (49) 

»,!!(£, 0 , 0 ) = -^£ 2 '-\ 2 N ,- Z 1 

v ' 1024ne L 7 12 

1 


- - ^4 hi b - hi a 


(53) 


^(v 2 + l),b=§(4v 2 
straightforward to check the following relation 


where a = |^(u 2 + 1), 6 = |^(4u 2 + 1). Finally it is 


2512'(0,0,0) = 


TTMe^m 8 
1024w 
— - ^4 In b — In a^j 




(54) 


While the probability of the pair production takes 
the form 


r 7T Me T m 8 

p = exp i ^02^r 

— ^ ^4 In b — In aj j 


2TV 7 — — 
7 12 

« 0. 


This end the proof of proposition (2). 


(55) 

□ 


IV. SOLUTION OF THE DIRAC EQUATION 


/Ci = i 7 0 do ~ *7 2 &2 — *7 3 fc 3 + i'f’Et + 9ioj^° 

—3257 1 — m(l — ujt), (56) 

IC 2 = — 7 1 <9i — i^ 2 Bx + mux. (57) 

The equation (12) takes the form (/Ci + £2)^(2, x) = 
0 and admit separate variables as ip(t,x ) = 

ij;(t)if)(x). For a constant A € C, we get the two 
eigenvalue equations 


ICii/j(t,x) = A i/)(t,x) 
K. 2 ip(t,x) = A ip(t,x). 


(58) 

(59) 


Consider the equation (58). We write the four 
vector ip(t) as ifi(t) = (ipi(t),ip 2 (t)) and i/>j(t) = 
(■ ipja(t ), ipjb(t)), j = 1, 2, and the equation (58) leads 
to 

(L'L + D 2 +CC')ip la (t) = 0, (60) 

(LL'+ D 2 + CC')^2a{t) = 0, (61) 

Cl/Jib{t ) = L'l/>2a(t) ~ Dlpi a (t), (62) 
Ci/j 2 b{t) = - Dlp2a(t ) (63) 


where 


Let 


L = i<9o + 9iw — m(l — uit) + A, 
L' = —id 0 — 9*w + to(1 — ut) + A 
C = —3w — k2, C' = —323 + Jt2, 
D = — ifc 3 . 


/(/) = 7V 2 + (~^( m “ A)w - 9w - fc 3 )t, (64) 


4 


<?(/) = (-DHiE)n 


mu 

2E^ m - X) 


(65) 


+ jp, (k 2 - (m - A) 2 + imuj - * (66) 

4 = {^k 3 ~ (m - A )muj. (67) 

The solution of the equations (60) are a linear com¬ 
bination of Hermite and (1,1)-hypergeometric poly¬ 
nomial given by 


4>ia(t) = c 1 e f{ - t) n\5 1 El 5 2 E +g{t) 


In this section we give the solution of the Dirac 
equation (12). We consider the operators K. 1 and 
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+ c 2 e f(t \F 1 l-^-,(S^ + g(t)) } 
i> 2 a{t) = c 1 e f( ' t) 'H <$£, 6% + g(t) 

+ c^.F^f^iSl + git)) 1 2 }. ( 68 ) 


the solutions of the equation (62) can be simple ob¬ 
tained using the followings identities: 

dH(a,b + ct) „ 

---= 2acH(— 1 + a, b + ct), (69) 

at 

d iFi(a, b , ct 2 + dt + e) 
dt 

_ a(d+2d) + + ^^2 _)_ fa + e)(70) 

o 

Now, consider the equation (59). Using the Dirac 
matrices (6), we get 


— (<7i<9i + ia 2 Bx)ip 2 {x) + (mixx — A)^i(x) = 0(71) 
(<7 i<9i -(- ia 2 Bx)ipi(x) + (mojx — X)ip 2 {x) = 0(72) 


where if)(x) = (ipi(x), ^(x)). Let us define the 
quantities b{B), r(B) and s(B) as 


b{B) 

r(B) 

s(B) 


_ 

2 B’ 

2mu)\(— 1) J 

V2(iB) i ’ 

= (-l)iv^(iB)^. 


For ipj{x) = ( ipja(x),ijjj b (x )), j = 1,2. The solu¬ 
tions of the equation (71) are 


4>ia(x) 


ipib(x) 


c{D 

c 2 D 


c{D 

c 2 V 


b{B),—r{B) + s(B)x 
b(B),—ir(B) + is(B)x 
1 — b(B), —r(B) + s(B)x 
1 — b(B), —ir{B) + is(B)x 


(73) 


(74) 


However, the equation (72) can be split into 

dl-B- B 2 x 2 + A 2 - 2mXuix tp la (x) = 0 (75) 
df + B — B 2 x 2 + A 2 — 2mXtjjx ipib{x) = 0, (76) 


and the solutions are well given by the following: 

ip 2 a(x) = -—^^{di + Bx)ipib, (77) 

A — mix x 

ip 2 b(x) = -——— (9i - Bx)ip la . (78) 
A — mixx 

where the identities 

—V(a,bx + c) = ~(bx + c)V(a,bx + c) 
dx 2 

— bT>(l + a,bx + c), (79) 


are usefull. 

V. CONCLUSION 


In this paper, we have computed the probability 
density of pair production of the fermion particles. 
The case of vanishing EM fields is scrutinized ex¬ 
plicitly. Hereafter we will shed light on the case of 
non-vanishing EM fields, which has not been entirely 
considered in this paper. In the other hand, we have 
solved the Dirac equation coupled with gravitation 
field, using the separation of variables. The solutions 
are expressed in terms of liypergeometric functions. 
The limit where the deformation parameter u> tends 
to zero is given. 
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